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^ , Abstract 

a decay and its relation with chiral phase transition are discussed at finite tem- 
perature and density in the framework of the Nambu-Jona-Lasinio model. The 
decay rate for the process a — > 2ir to first order in a 1/N C expansion is calculated 
. as a function of temperature T and baryon density rib- In particular, only when the 

!> \ chiral phase transition happens around the tricritical point, the a decay results in 

a non-thermal enhancement of pions in the final state distributions in relativistic 
. heavy ion collisions. 

00 ; 

It is generally believed that there are two QCD phase transitions in hot and dense 
nuclear matter []I||. One of them is related to the deconfmement process in moving from a 
hadron gas to a quark-gluon plasma, and the other one describes the transition from the 
^ ■ chiral symmetry breaking phase to the phase in which it is restored. 

At low energies, there is no obvious reason for the presence of the scalar meson a in the 
study of chiral properties. However, there is a definite need0 for a at finite temperature 
and density in the investigation of chiral symmetry restoration. For instance^, the 
| contribution from a to thermodynamics can for most case be neglected, because of its 

heavy mass, while it plays an important role in the region around the critical point of 
the chiral phase transition, since there, m a = 0. Because it is still difficult to extract 
definite information from the lattice simulations with nonzero baryon density, we need 
QCD models to investigate the phase transitions at finite density. From the study of 
several models Q, the chiral phase transition in the chiral limit with massless pions is of 
second order at high temperature T, and of first order at high density rib. Therefore, there 
is a tricritical point P in the T —rib plane which separates the first- and second-order phase 
transitions. In the real world with nonzero pion mass, the second-order phase transition 
becomes a smooth crossover and the tricritical point P becomes a critical end-point E of 
a first-order phase transition. In relativistic heavy ion collisions, for the choice of control 
parameters like colliding energy and impact parameter such that the freeze-out of the 
system occurs near the point E, a is one of the most numerous species at the freeze-out, 
since it is nearly massless. During the expansion of the system, the in-medium a mass 
rises towards its vacuum value and eventually exceeds the im threshold. As the ann 
coupling is large, the decay proceeds rapidly. Since this process occurs after freeze-out, 
the pions generated by it do not have a chance to thermalize. Thus, it is expected that 
the resulting pion spectrum should have a non-thermal enhancement at low transverse 
momentum. 
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In this letter we discuss the decay process a — > 2tt at finite temperature and density in 
the framework of the Nambu-Jona-lasinio (NJL) model which is believed to describe the 
chiral properties well. This was firstly studied at zero density by Hatsuda and Kunihiro|Q. 
Our motivation here is to see that in which temperature and density region of freeze-out 
there is a non-thermal enhancement of pions due to the a decay, by calculating the decay 
rate in the whole T — nj, plane. 

The two-flavor version of the NJL model is defined through the lagrangian density, 

L NJL = - m )^ + G[(^) 2 + $iW) 2 ], (1) 

where only the scalar and pseudoscalar interactions corresponding to o and it mesons, 
respectively, are considered, ip and ^ are the quark fields, r is the SU(2) isospin generator, 
G is the coupling constant with dimension GeV~ 2 , and mo is the current quark mass. 

In an expansion in the inverse number of colors, 1/N C , the zeroth order (Hartree) 
approximation for quarks together with the first order (RPA) approximation for mesons 
gives a self-consistent treatment of the quark-meson plasma in the NJL model. To the 
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lowest order, namely (1/N C ) , the Feynman diagram for the decay process a — > 2tt is 
sketched in Fig.l. The partial decay rate in the rest frame of a is given in terms of the 
Lorentz-invariant matrix element M by 

%^ = ^4W> ( 2 ) 
ail 32tt z mi 



where |p| = y^j- — vn\ is the pion momentum, m a and m n are respectively a and it 
masses, using an obvious notation for the M-matrix element, one has 



M = g a g 2 A a7T7T (3) 



with 



A™ = TJ2 e lUnr > f ? ^%Tr5(q, iw n )T v S(,q + p, iu n + ^)r^(q, iu n + m a ) . (4) 

J (ITT) 6 2 

Here u n = (2n + 1)ttT, n = 0, ±1, ±2, ... are fermionic Matsubara frequencies, = i^y^r 
is the pseudoscalar interaction vertex. The quark propagator is denoted by S(q, iuj n ) = 
(7o(ia; n ) — 7 • q + m q ) / (iiuo n ) 2 — E 2 ^, with E 2 = m 2 + q 2 . In Eq.(||), Tr refers to the 
trace over color, flavor and spinor indices. The pion-quark and sigma-quark coupling 
strengths in Eq.(|3]) are determined in the model via|| 



f71I ,(/, (1 .0: V. //) 



dkl 

dU a {k ,0;T, fJl ) l 



\kl=m\ 1 



9a (T,fl) = ^2 ' |fcg =m 2 , (5) 



where II,,- and U a are the standard mesonic polarization functions for tt and a. The 
dynamically generated quark mass m q (T, fi) and the meson masses m n (T, /1) and m a (T, /i) 
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are calculated using the usual gap equations || in the Hartree approximation for quarks 
and RPA approximation for mesons. After evaluation of the Matsubara sum in Eq.(|4]), 
one has in 



8(q • p) 2 - (2m 2 + 4m 2 )q • p + m 2 /2 - 2m 2 £ 2 



ml — AE 2 ) ((ml — 2q • p) 2 — m 2 i? 2 



(6) 
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where /f(x) = (l + e x ^ T J is the Fermi-Dirac distribution function. 

a can decay into neutral and charged pions. By considering the exchange contribution 
and < 9 < n /2 for a — > 2ttq due to the identical particle effect, the total decay rate can 
be written as 

F a ->2ir (T, /i) = r o -_»27r (T, ^) + r o -_„ r+7r _ (T, /i) 



3X/7 m ^l^(T,/.)| 2 (l + M^)) Z , (7) 



?7r ml 



where we have taken into account the Bose-Einstein statistics in terms of the distribution 
function fs(x) = (e x ^ T — lj for the final state pions. 

From the comparison of a— it— 7r coupling strength g CT7r7r defined through the lagrangian 
Li ~ ^o-tttt " 71 " ' 7T and the M-matrix element (^|) in the NJL model, we have 

<W(T, //) = 2M(T, /i) = 2(^A(T, //) . (8) 

Now let us turn to numerical calculations of the coupling strength g a7Tn and the decay 
rate r o -__ >2 7r as functions of temperature T and baryon density n b . The relation between 
rib and the chemical potential /i is rather complicated and is also a function of T, 



N c N f f d 3 q 



(27T) 



tanh — (E g + n) - tanh ^ - //) 



(9) 



where the factor 3 reflects the fact that 3 quarks make a baryon. In the calculation we 
use the dynamical quark mass m q = 0.32 GeV, the pion decay constant /„- = 0.093 GeV 
and the pion mass m n = 0.134 GeV in the vacuum (T = ji — 0) as input to determine 
the 3 parameters in the NJL model, namely the coupling constant G = 4.93 GeV~ 2 , the 
momentum cutoff A = 0.653 GeV and the current quark mass mo = 0.005 GeV. 

We first discuss the temperature dependence of g an7T and r cr ^ 2 7r in the case of zero 
baryon density which corresponds to the central region of relativistic heavy ion collisions. 
The coupling strength g an7T is about 2 GeV in the vacuum and almost a constant in 
the temperature region T < 0.1 GeV, and then decreases smoothly with increase of the 
temperature, as shown in Fig. 2. The energy condition for the decay process o — + 2n to 
happen is m G > 2m n in the rest frame of a. At the threshold temperature T ~ 0.19 GeV 
determined by m a (T) = 2m 7r (T), g aiT1T has a sudden jump from about 1.2 GeV to zero. 
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In a wide temperature region, the decay rate (0) goes up gradually from its vacuum 
value ~ 0.09 GeV to the maximum ~ 0.11 GeV at T ~ 0.17 GeV, and then it drops 
down rapidlly but continuously to zero at the threshold temperature. The fact that the 
maximum value of T is not in the vacuum, but at T ~ 0.17 GeV is beyond our expectation. 
It originates from the important Bose-Einstein statistics of the final state pions at high 
temperature. This can be seen clearly in Fig. 2 where we show the decay rates with and 
without consideration of the Bose-Einstein statistics. The rapid decrease of T a ^ 2 n at high 
temperature is due to the behavior of g an7T and the phase-space factor Jm^/A — m^jrv? a . 

The density effect on g an7T and r cr _ +2 7r is shown in Fig. 3. Different from the temperature 
effect in Fig. 2 where the statistics of the final state pions is a crucial factor, there is no 
difference between with and without such statistics when we discuss density effect at low 
temperature (see Fig. 2) which corresponds to the region of the first-order chiral phase 
transition in the T — plane. Instead of the statistics, the chiral restoration dominates 
the behavior of g U7Tn and r cr _ +2 7r- At T = 0, the mixing phase of the first-order chiral 
transition is in the density region 0.33 < n^/no < 1.56 where no is the baryon density in 
normal nuclear matter. At the up limit rib/no = 1.56, m a < 2m n , there is no a decay, 
while at the low limit nb/n = 0.33, m a > 2m v , a decay begins with maximum rate. 
With the decrease of the density, the decay rate and the coupling strength drop down 
from the maximums to their vacuum values, respectively. The fact that the maximum 
decay rate and the maximum coupling strength local at the low density limit of the mixing 
phase is directly due to the first-order chiral transition. The situation for T = 0.02 GeV 
is similar, but the mixing region becomes narrower, 0.59 < nf,/n < 1.44, and the jump of 
the decay rate (coupling strength) becomes smaller. As the temperature increases further, 
the density interval of the mixing phase becomes narrower and narrower, and the jump 
becomes smaller and smaller. Finally when approaching the end-point E of the first-order 
transition at T ~ 0.03 GeV, the interval and the jump disappear. The decay begins at 
the point E and proceeds faster and faster with the decrease of the density. The decay 
rate goes up from zero at the phase transition point E to the maximum at low density. 
While the coupling strength has still a jump at the point E, the decay rate is continuous. 

To see the parameter dependence of the above discussion, we considered another set 
of parameters, G = 5.46 GeV~ 2 , A = 0.632 GeV and m = 0.0055 GeV. It makes shifts 
of the magnitudes of g a7T7T and T a ^ 2 n, but does not change their shapes, especially the 
behavior around the threshold temperature or the chiral phase transition point. 

From the above investigation, we see that when a decay into 27r happens in low density 
or low temperature region, the maximum decay rate does not locate in the vacuum, but 
is very close to the threshold point or just at the first-order chiral transition point. Only 
when the decay trajectory passes through the region around the end-point E of the phase 
transition, the position of the maximum decay rate is far from the point E where the 
sigma has minimum mass and therefore can be numerously produced. According to this 
picture, when the system of a relativistic heavy ion collision evolves in the low density or 
low temperature region, the a decay begins with maximum rate immediately after or even 
simultaneously at its starting point of creation. Since the a decay happens already at 
the freeze-out, there is no obvious non-thermal enhancement of pions. However, when the 
freeze-out occurs around the tricritical point, the a decay happens far from the freeze-out, 
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there is really a non-thermal enhancement of pions which can be considered as a direct 

signature of chiral phase transition. 
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Figure Captions 



Fig.l: The Feynman diagram for a decay into 2n to the lowest order in 1/N C expan- 
sion. The solid lines denote quarks and antiquarks, the dashed lines denote mesons. 

Fig. 2: The total decay rate r cr _ >2 7r and the coupling strength g a7T7r as functions of the 
temperature at = 0. The thick and thin solid lines indicate the decay rates with and 
without consideration of Bose-Einstein statistics for the final state pions. 

Fig. 3: The total decay rate r cr __ +2 7r (solid lines) and the coupling strength g a7T7T (dashed 
lines) as functions of the baryon density scaled by no for three values of the temperature. 
The black circle and black lozenge indicate the starting- and end-point of the mixing 
phase at T = and T = 0.02 GeV, respectively. 
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